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The system-plus-reservoir (SPR) model is the most common and effective approach to study 
quantum dissipative effects. Indeed, it makes quantization possible by considering the whole energy- 
conserving system, while the reservoir's degrees of freedom, assumed to be harmonic, can be traced 
out by the path-integral technique, leading to a formulation that only includes the system of interest. 
In the standard SPR model the environment is only coupled with the system's coordinate and turns 
out to quench its quantum fluctuations. However, there are physical systems coupled with an 
environment whose 'coordinates' and 'momenta' can be completely interchangeable (e.g., magnets), 
so an SPR coupling must symmetrically affect both canonical variables. In this paper such a general 
environmental coupling is studied in the case of a harmonic oscillator. It is found that quantum 
fluctuations are generally enhanced by environmental coupling, with an unexpected nonmonotonic 
behavior. This leads one to speculate about the possibility that spin-lattice coupling could drive 
the 2D Heisenberg antiferromagnet to reach its quantum-critical point. 



I. INTRODUCTION 

The most popular approach to the theory of quantum 
dissipative systems [l| involves the system-plus-reservoir 
(SPR) model where dissipation is assumed to arise 

from the interaction with a reservoir (or environment, 
or bath); the latter is constituted by very many degrees 
of freedom which are individually weakly perturbed by 
the system, so that (a) the reservoir can be assumed to 
remain at thermal equilibrium and (b) its degrees of free- 
dom can be modeled as harmonic oscillators, which are 
(c) linearly coupled with some system's observable. This 
observable is usually taken as a function of the system's 
coordinate only, and the case of strictly linear dissipation 
corresponds to the case of a linear function. 

Within such a model it is possible to show that the (re- 
tarded) dynamics of the system, once the bath variables 
arc eliminated, is indeed described by a Langevin equa- 
tion Q, where the dissipative memory function is con- 
nected to the characteristics of the bath and its coupling 
to the system. This contact between the microscopic de- 
tails of the bath and the phenomenological description is 
valuable, as it allows to connect with the phenomenol- 
ogy the effects onto the quantum thermodynamics of the 
system. 

Originally, the picture of a more or less fictitious har- 
monic bath was meant to give a physically sound tool for 
describing quantum dissipation: a Hamiltonian quantiza- 
tion is made possible because the overall system includ- 
ing the bath is isolated. However, in several real cases it 
is possible to precisely devise the physical environment 
and its microscopic Hamiltonian. For instance, a mag- 



netic system is intrinsically coupled with its underlying 
lattice; this is demonstrated by the observed occurrence 
of effects that can be as dramatic as to yield a lattice 
distortion (spin-Peierls transition 0, Q). Now, while 
the motion of the lattice ions can be fairly described in 
terms of linear excitations (phonons), it is also possible 
to microscopically model the coupling between the spins 
and the ion positions, so that the vibrating lattice can 
be eventually treated as an environmental bath for the 
spins. Within this picture one could expect and theoreti- 
cally calculate observable modifications of the behavior of 
quantum magnets, such as a shift of the predicted critical 
temperature |9l4ll| ■ As spin Hamiltonians are often dealt 
with through a spin-boson transformation reducing spin 
operators to ordinary canonical variables, thinking of a 
possible environmental coupling would symmetrically in- 
volve both the coordinates and the momenta, with no 
privileged role: this has been shown in Ref. [12| for the 
case of the easy-axis XXZ magnet. Therefore it is neces- 
sary to generalize the concept of the dissipative system 
in a much more general way than it was done when con- 
sidering anomalous dissipation [l^Hl5|. where the bath 
is assumed to be coupled with the momentum. 

Such a generalization has been used for studying the 
coherence properties of a magnetic impurity in a ferro- 
magnetic environment : it has been shown that some 
quantum effects of the two baths partially cancel, leading 
to a persistence of coherence that has been called quan- 
tum frustration of decoherence [l7|, and this motivated 
for a more detailed study of the dynamics of the quantum 
oscillator coupled with two baths |lq |. 

The aim of the present paper is to provide a deeper 
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understanding of the thermodynamic behavior that fol- 
lows from the intrinsic dynamical character of quantum 
thermodynamics. Indeed, besides the dynamical features 
studied in Ref. [18], one expects that also the system's 
static quantities could exhibit effects that witness for 
quantum frustration. The present work gives an overview 
of the different - sometimes unexpected - features that 
a generalized environmental coupling can produce. 



From previous studies it is well known that, roughly 
speaking, the quantum fluctuations of the variable to 
which the bath is attached are quenched and those of 
its conjugate variable are enhanced, in such a way that 
the Heisenberg uncertainty principle holds. However, it 
is not clear which would be the prevailing effect when at- 
taching one or two independent baths to both canonical 
variables, especially when the coupling has an equiva- 
lent weight onto coordinate and momentum (symmetric 
coupling). We are going here to answer such questions 
by assuming the simultaneous existence of environmen- 
tal coupling with both canonical variables, taking an ex- 
actly solvable reference system, i.e., the harmonic oscil- 
lator. We will find that the presence of a symmetric 
environmental coupling generally enhances the quantum 
fluctuations of coordinate and momentum, but that this 
does not necessarily occur in a monotonic way with the 
coupling strength. Moreover, as such phenomena can 
show up in a qualitatively similar way in nonlinear sys- 
tems, one can think of the possibility of finding reentrant 
behavior driven by the environmental coupling strength, 
e.g., when the system is close to a (quantum) phase tran- 
sition. 



In Section |TT] the framework of general environmen- 
tal coupling is introduced considering two independent 
reservoirs coupled to the coordinate and to the momen- 
tum, while in Section [TTT] the alternative case of one single 
bath coupled to both canonical variables is approached; 
there, the main results for the environmental effects upon 
the mean-square fluctuations of the harmonic oscillator 
are reported and discussed for coupling with the only 
environment coordinates and also when the environmen- 
tal momenta are involved. In Section IIVI we summa- 
rize the results and draw some conclusions and specu- 
lations about the possible implications of symmetric en- 
vironmental coupling. In Appendix [A] the backbone of 
the standard theory of quantum dissipation is briefly re- 
viewed by recalling the definition of relevant quantities 
and setting the formalism adopted in this paper; Ap- 
pendix [B] reports the details of the evaluation of the gen- 
eral Gaussian path integral repeatedly used throughout 
the paper. 



II. COORDINATE AND MOMENTUM 
COUPLED WITH TWO INDEPENDENT 
ENVIRONMENTS 

We consider as system of interest a quantum harmonic 
oscillator, 

H s = \(a 2 p 2 +b 2 q 2 ), (1) 

where we prefer to introduce two parameters a and b 
in the place of the commonly used mass m = a~ 2 and 
frequency ujq — a b, in such a way to emphasize the sym- 
metry between p and q. In addition, we assume h=l, 
i.e., the commutator [q,p] = 1. 

A. Influence action and fluctuations 

As motivated in the Introduction, we wish here con- 
sider the case of two baths, one coupled with the coor- 
dinate and one with the momentum. The known results 
summarized in Appendix [A] tell us that the two indepen- 
dent baths are in conflict, so the first natural but nontriv- 
ial question is which effect prevails onto the mean-square 
fluctuations of, say, the coordinate. Of course, the answer 
also depends on the relative intensities of the couplings; 
however, if the bath couplings are identical, the behavior 
is not easily foreseeable. 

We start then by coupling the canonical variables of the 
oscillator Hamiltonian (flj with two independent baths, 
i.e., 

Q 

+ \Y,[4(h-p) 2 + bl<iV\ ; (2) 

here the subscripts are used to distinguish between the 
two independent sets of bath parameters {a a ,b a } and 
{ap, bp}, corresponding to the frequencies {ui a — a a b a } 
and {ujp = ap bp}. In analogy with Eq. (|A4I) we are thus 
lead to introduce two 'memory functions', J p (t) and 7 g (t), 
whose Laplace transforms read 

Or b 2 

p z +^p « * + 

The equations of motion for the system's coordinate 
can be again put into the form of a Langevin equation, 

'q\t) + J dt' T(t') q{t-t') + LO 2 q(t) = F(t) , (4) 

where 

r(t) = b 2 lp (t) + a 2 lq (t)+j pq (t) (5) 

and 

7 M (*)= / dt ' 7 P (07,(f-0 ■ ( 6 ) 
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Assuming j q (t) and 7p (t) to vanish for t < and to be 
positive decreasing functions for t > 0, then T(t) also 
satisfies the condition of causality, but as j pq (t) < the 
positivity of T(t) is not guaranteed. Basically, this means 
that there isn't a standard phenomenological dissipa- 
tive counterpart of our system-plus-two-reservoirs model. 
Therefore, it is better to refer to the two-bath coupling 
as environmental, rather than dissipative, coupling. Of 
course, in physical applications the baths and their inter- 
action with the system have to be microscopically charac- 
terized. Anyhow, since in the limit where one of the baths 
can be disregarded we expect a (standard or anomalous) 
dissipative behavior, we argue that a phenomenological 
form for both memory functions is physically reasonable; 
we will make use of it later on. 

The calculation of the influence action can be made 
separately for the two baths in the very same way that 
leads to Eq. (|X2"]h 



S\[p,q] 



1 

2 jo 

P 



dudv! \k p {u—u')p(u)p(u') 

+k q (u-u') q(u)q(u')] 



fluctuations 

VP /cl = Jtf. = J ' W )cl = ^2 = £^2 S ( H ) 

noting that the classical contribution is unaffected by the 
environment, it is useful to separate the remaining purely 
quantum contribution, which includes the whole environ- 
mental effect, 

/ ~2\ _ ^ ^-v fo 2 + k qn 

'™ ~ J ^ vl + (a 2 +k pn )(b 2 +k qn ) ' 



<« a > pq = 5E 



(12) 



^ A ^jJTt PnP—n ~t~ C[nQt — 7i) 



(J) 



For the isolated oscillator these quantities reduce to 

x ^ x Pq 2a V 2 /?w / 

<0») = ^( coth ^£_M (13) 
vy ; pq 26 V 2 /3wo/ 

which look more familiar by noting that b/a — muiQ. 



where the Matsubara components of the kernels k p (u) 
and k q (u) read 



kpn — ^ 



2 2 



K 9» Z^„2,,,,: 



(8) 



so that they are connected with the memory functions 
7 P (z) and J q (z) by a relationship analogous to Eq. (|A5|) . 

Eventually, including the isolated system's action and 
using the general result of Appendix [B] one finds the par- 
tition function 



/3W0 ^ ^ + (a 2 + fcpn)(& 2 + fc 9 n) 

and the mean-square fluctuations 

(f) = tX 



(9) 



& 2 + k qn 



p ^ + (a 2 +k pn )(b 2 +k qn ) ' 



'■pit 



(10) 



which, as expected, are symmetrically related and reduce 
to the known forms in the limits of standard (k pn — > 0) 
and anomalous (k qn — > 0) dissipation. The summations 
in Eqs. (|10[) are convergent, provided that for large n both 
kernels are such that k n /n —> 0. Together with the result 
(pq + qp) = 0, Eqs. (fTUj) fully determine the Gaussian 
(reduced) density matrix corresponding to the oscillator 
Hamiltonian ([1} plus the environmental interaction @. 

The zeroth Matsubara component of Eqs. (fTU)) coin- 
cides with the classical contribution to the mean-square 



B. Environmental effects 

Let us consider the general expressions (|12l) for the 
purely quantum mean-square fluctuations. Due to the 
interchangeability of the canonical variables, it is suffi- 
cient to study, say, the coordinate fluctuations only. 

To proceed, we take the minimal phenomenological 
form for the memory function of each bath, i.e., the 
Drude model [J, defined by 7 (i) = 77-- 1 e~ t/r o and 
characterized by the intensity 7 and the memory time 
r D . Its Laplace transform, to which the bath's kernel is 



(14) 



related via Eq. (|A5j) . reads 

7(2) 



The Drude cutoff frequency oj d ^t^ 1 characterizes the 
environment and its interaction with the system, so it 
is expected to be of the order of the Debye frequency 
in the case of a phonon bath. In the limit of vanish- 
ing memory time, i.e., large w D , one recovers the so- 
called Ohmic (or Markovian) model, which is memory- 
less, j(t) = 7<5(t— + ). Actually, a finite Drude frequency 
is necessary in order to get finite results for both quan- 
tities in Eqs. (fl2l ; on the other hand, it is indeed known 
that Ohmic response is incompatible with quantum me- 
chanics (I9ll20j. 

The Drude kernels we are going to use are then 



^ P n 7p 



1+T D K| 



k, 



qn 7g 



1 + T D K| 



(15) 



where for both kernels the same memory time r D is cho- 
sen for simplicity, while the amplitudes 7 P and j q will be 
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varied independently. In addition, without loss of gen- 
erality we can take a symmetric form for the oscillator 
Hamiltonian ([T]), i.e., a 2 — b 2 — ljq; in this way both 
7 g and 7 P are dimensionless. Taking lo$ as the frequency 
unit leads to a natural dimensionless formulation, with 
the reduced temperature t — l/(/?wo), the Matsubara fre- 
quency u n = 27rfnu;o, and so on; all this is tantamount to 
setting too = 1, and in such units the classical part of the 
dimensionless fluctuations is equal to the temperature, 

Therefore, the mean-square fluctuation we have to 
study is 



(<f ) pq = £*(*,7p,7?,t"d) , 
which has the dimensionless explicit expression 



(16) 




FIG. 1: Mean-square fluctuation of the coordinate at t = as 
a function of both 7 9 and 7 P , for r D = 0.01. 



2 *E 



1 + k n 



< x {2ntn) 2 + (1 + k pn )(l + k qn ) ' 



with 

kpn 7p 



2ntn 



1 + T n 27rfn 



kqn — 7q 



2irtn 



1 + Tx, 2irtn 



(17) 



(18) 



Note that a nonzero memory time r D is required to make 
the series (fTT)) convergent. 



1. Zero temperature 

Taking the limit t — > 0, one sets x — 2ittn and dx — 2nt 
to write Eq. (fT7|) as an integral, 



where 



dx 



,, 7T a; 2 + [l+7 p g(x)][l+7 g g(a;)] 

(19) 




FIG. 2: Mean-square fluctuation of the coordinate at t = as a 
function of 7 9 for fixed values of y p . As in the case of standard 
dissipation, 7 P = 0, the effect of increasing the intensity of the 
coordinate coupling is a decrease of (q 2 )- 



1 + T^X 



(20) 



From Eq. (jT9")) one can see that the result of the isolated- 
system agrees with Eq. (|13p in dimensionless form, 



a(0,0,0,r D ) = 



1 



dx 



c 2 + 1 2 ' 



(21) 



and also recover a known analytic result [T| in the case 
of Ohmic standard dissipation, 



a(0, 0,7,0) 



2 tan-VlEZ (7<2) 



7T-\/4-7 2 



7 



2 tanh" 1 ^tLA ( 7 >2). (22) 



7T\/7 2 -4 
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In the general case Eq. ([19]) has to be evaluated numer- 
ically. This is what we have done and Fig. [1] shows the 
three-dimensional resulting plot for a(0, 7p , 7g , r D =0.01); 
its shape gives an overall idea of the combined effect of 



two baths. However, in order to highlight some impor- 
tant features which are barely apparent, it is necessary 
to analyze its cross-sections. 

In Fig. [2] we first consider how a varies while increasing 
the coordinate coupling j q : also for nonzero 7 P the effect 
is a decrease of a, as if the g-bath were a position measur- 
ing device, causing localization of the particle. From the 
same figure, and more clearly from Fig. [3l the comple- 
mentary interpretation can be given to the action of the 
p-bath, generalizing the observations made in the case 
of anomalous dissipation (l4l. Il5|; note, however, that as 
soon as 7 9 is switched on the steepness of the rise of a is 
considerably suppressed. 

Given the opposite effect of the two baths, it is natu- 
ral to ask what the result will be when both have equal 
strength, j p = j q : this is a nontrivial issue, especially 
from the point of view of the application to real prob- 
lems where the canonical variables play symmetric roles, 
as mentioned in the Introduction. 

The answer is given by the curves reported in Fig. 2) 
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FIG. 3: Mean-square fluctuation of the coordinate at t — 
as a function of 7 P for fixed values of j q . As in the case 
of anomalous dissipation, 7 9 = 0, the effect of increasing the 
intensity of the momentum coupling is an increase of (q 2 )- 
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FIG. 4: Mean-square fluctuation of the coordinate at t = as 
a function of 7 = 7 g = y p , i.e., for equally acting environments, 
for selected values of the Drude time r D . It appears that be- 
tween the competing effects of the two baths, the prevailing 
one is that of increasing (q 2 ); however, as the coupling inten- 
sity increases a maximum is reached and followed by a slow 
decrease towards the isolated-oscillator value. 

a = (<? 2 ) pq is larger compared to the isolated oscillator, 
i.e., the quantum fluctuations of both canonical variables 
are enhanced by a symmetric environmental coupling. 
However, an unexpected behavior shows up while rising 
the intensity of the environmental coupling: after an ini- 
tial rise the mean-square fluctuation shows a maximum, 
whose location and intensity depends upon the Drude 
memory time, and eventually slowly decreases back to- 
wards the isolated oscillator value a = 1/2. 

It is obvious that, as the size of the zero-temperature 
fluctuations is affected by the mechanism of environmen- 
tal coupling, one could formally modulate its intensity 
7 to drive through a quantum phase transition (QPT) a 
system that is close to it. At difference with standard 
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FIG. 5: Mean-square fluctuation of the coordinate for r D = 0.1 
(a) and t d =0.01 (b) at t = along the line "f P ~cy q as a 
function of j q , for different values of c. The dashed line, for 
c= 0.2214 in (a) and c= 0.1193 in (b), separates the decreasing 
from the nonmonotonic behaviors. 



dissipative systems, e.g., resistively shunted Josephson 
junction arrays, where a stronger dissipation leads to- 
wards the 'ordered' (superconducting) phase [2l|, in the 
generalized case a higher 7 would most commonly in- 
crease disorder. One can also speculate about the fact 
that the nonmonotonicity displayed in Fig. 2] opens the 
possibility to observe reentrant behavior around a QPT, 
namely by further rising 7 beyond the critical value one 
could observe the restoration of the isolated-system's or- 
dered phase. 

Having ascertained that the combined effect of the two 
baths is opposite to that of the single q-bath, it is nat- 
ural to look for the evolution one observes by gradually 
switching the p-bath coupling on. In Figs.[5]the behavior 
of a is followed along straight lines in the plane (79, 7p), 
i.e., for fixed ratios c = 7 P /7 g between the intensities of 
the two environmental couplings. These figures, which 
differ for the choice of the Drude time, are essentially 
similar and display the competing interplay between the 
intensities of the environmental couplings for q and p: if 
the latter is small enough, the quenching effect of the for- 
mer prevails. Increasing 7 P for any fixed 7 9 , a must cross 
the isolated-system value a = 1/2 for some 7 P : along the 
line 7 P (7 9 ) the competing effects cancel each other and 
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FIG. 6: The function 7 P (7q) corresponding to the cancela- 
tion of the competing effects of the two baths, so that the 
mean-square fluctuation of the coordinate keeps the isolated 
system's value. The curves refer to t~0, so (q 2 ) = 1/2, and 
to selected values of the Drude time r„ . 



the mean-square fluctuation of the coordinate is station- 
ary; since a > 1/2 for 7 P = j q it follows that 7 P (7 g ) < 7 9 - 
From Eq. (flU| one can analytically obtain the initial slope 
of these curves as a function of the Drude memory time 



%(0) 



2-irr n -4r2lnr n + 2t 2 



■ 7TT~ 



-41nr n - 2 + 37rr n - 2t 2 



(23) 



which turns out to vanish in the Ohmic limit and tends 
to 1 for large r D . This is in agreement with the numerical 
calculation of 7 P (7<j) reported in Fig. El 

Fig. [5] also suggests that, keeping the ratio c = 7 P /7 9 
fixed, for large coupling intensity (q 2 ) tends to a finite 
value; one can indeed find that 



(24) 



Substantially, it turns out that, after a possible initial 
increase of the spread of one of the canonical variables 
when the interaction is switched on, a very strong en- 
vironmental coupling eventually tends to localize again 
both p and q, as the above limit also entails that 




(25) 



the shrinking of the fluctuations is visibly more effective 
for the variable with larger bath-coupling strength, but 
the product of the fluctuations obeys the constraint of 
the uncertainty principle, 



(26) 



In other words, a very large environmental coupling leads 
to a situation where the uncertainty is squeezed along one 
axis in the p-q plane with respect to the symmetric spread 
in the free system. 



FIG. 7: Pure-quantum mean-square fluctuation of the coor- 
dinate for t d = 0.1 and different values of 7 = y q — j p (i.e., for 
equally acting environments) as a function of t. 
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FIG. 8: Pure-quantum mean-square fluctuation of the coor- 
dinate at different temperatures as a function of 7 = 7<j = 7p , 
for t d = 0.1. The dashed curves for t > 2 represent the large-f 
analytic approximation (|27|) valid for t ^> 1/(2ttt d ) ~ 1.6. 



2. Finite temperature 

As Eq. (fTT)) shows, the effect of environmental cou- 
pling only affects the purely quantum part of the fluctu- 
ations, which decreases with increasing temperature: so, 
one expects a softening of the effects observed at zero-T 
when T rises. However, the fate of the reentrant fluctua- 
tion amplitude displayed in Figs. 2] and [5] is not obvious. 
To unveil such issue we plot in Fig. [7] the behavior of 
a = (q 2 ) P q vs temperature for the situation of equally 
acting environments. It appears that a stronger envi- 
ronmental coupling weakens the decrease of (q 2 ) P q with 
temperature, and that for t > 1 a monotonic behavior is 
apparently restored. This fact is more evident by look- 
ing at the curves for increasing temperatures reported in 
Fig.[U 
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When the temperature is larger than the Drude cutoff 
frequency, t ^> 1/(27tt d ), the expressions (fT5)l for the ker- 
nels can be well approximated by their asymptotic values 
k qn ~ lq/T D and k pn ~ 7p/t d , so their dependence on n 
disappears and the summation (|17|) becomes analytic, 



(27) 



with Cp ? = 1 + 7p,g/T D : this approximation is shown 
as dashed lines in Fig. [5] for c p c q /2t^> 1, i.e., under the 
combined condition y / 7^7q 3> 2t d i > 1 /7r, the value of a 
tends to the constant c p /2c q . As the validity of this ap- 
proximation entails the disappearance of the maximum, 
the nonmonotonic behavior is expected for t < 1/(2itt d ). 



III. COORDINATE AND MOMENTUM 
COUPLED WITH A SINGLE ENVIRONMENT 

It can also happen that there is one single environment 
coupled with both canonical variables: the most general 
influence Hamiltonian can be found in Ref. [llj]. How- 
ever, there are two main scenarios of physical relevance: 
the first occurs when both q and p are coupled with the 
same bath variables (say, the coordinates), while in the 
second the q a 's are coupled with q and the p Q 's with 
p. The first scenario has a realization in the coupling 
of spins with lattice vibrations due to their modulation 
of the exchange integral; in the second case a natural 
model could be a spin-orbit interaction where the orbital 
angular momenta are thermalized by lattice vibrations. 



A. System coupled with bath coordinates 

Starting from the coupling Hamiltonian, 

& = lJ2 KpI + - c ^ - d «p) 2 \ > ( 28 ) 



we get the influence action: 

S\[p,q] = -- ^ [k qn q n q-n + k P n PnP-n + 2K n p n q-n] , 
n 

(29) 



with the kernels 

h 2 r 2 7/ 2 h 2 A 2 7/2 

, \ " ° a c a v n , \ "* ° a a a v n /on\ 

n 1 a - rj 1 



n ' a 



b 2 c d v 2 

"q L a u a^7i 



(31) 



It is immediate to check that the last kernel vanishes if 
c a d a = for any a, yielding the previous Eq. ([7]): in- 
deed, if the coefficients c a and d a are not simultaneously 
nonzero, the bath can be split into two independent baths 



and the two-bath case studied in the previous sections is 
recovered. 

In the case when all the coefficients c a and d a are pos- 
itive the spectral density [l[ J pq (uj) for k„, such that 



2vt 



J p q{w) 



(32) 



can be expressed in terms of those for the first two ker- 
nels, 

a 

JpH = ^bidiuJaS^-ujl) , (33) 

a 

as their geometric average as their geometric average [HI 

(34) 



Jpq{i0) = J J q (uj) J P (U)) . 



The baths are independent if the supports of J q (ijj) and 
Jr>(u)) have no intersection, i.e., if J pq (uj) vanishes for any 

Eventually, using the result of Appendix [B] the ex- 
pressions for the fluctuations in the case of the bath cou- 
pling (f28)) are found, 



oo 

^ n 



1 n=l 

2 °° 



- 2 




- (a 2 +fc pn )( &2 + 


kqn ) 


- K 2 : 










tv 2 - 


f (a 2 +fcp„)(&2- 















(35) 



(36) 



P n=l V n + (a 2 +k pn )(b 2 +k qn ) - nl • 



(w+#) p =-^E 



^1 + (« 2 +V)( &2 + V) - < ' 

As Eqs. (I3U1) and (|31l) imply the inequality k qn k pn > 
k 2 , it follows that the denominators are always positive. 
Note that if the bath is coupled to the variable p ± q 
(i.e., c Q = ±d Q ), then ±k„ = fc pn = fc g „ > 0: from the 
above variances it follows then that piq becomes smaller 
(larger) when the environmental coupling is switched on, 



<(P±<?) 2 > 



2 \ a 2 + b 2j r{k pn + k qn =p 2n n ) 



Pq ^ < + (a 2 +k P nW+k qn ) - «2 ' 

(37) 

which corresponds to a squeezed Gaussian distribution 
along a diagonal in the p-q plane, in agreement with the 
Heisenberg uncertainty relation. This observation also 
confirms the overall picture that the environment acts 
as a measuring device for whatever variable it is coupled 
with. 

Let us consider now a phenomenological guess fo the 
first two kernels (|3U|) . while the third one, K n , follows 
from Eqs. (|32p - (|3"4")) . Taking Drude kernels with equal 
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FIG. 9: System coupled with the bath coordinates. Mean- 
square fluctuation of the coordinate, Eq. (|40[) . at t = along 
the line "/ p =cy q as a function of j q , for different values of c. 
The dashed lines correspond to the curves for the two-bath 
case of Fig. [5] (b). 



B. System coupled with bath coordinates and 
momenta 

In this case the most general interaction Hamiltonian 
would take the form 



(42) 



note that in this case the coefficients c a cannot be ab- 
sorbed by means of linear canonical transformations of 
the bath variables. 

The calculation of the influence action leads to an ex- 
pression similar to Eq. , but for the appearance of an 
antisymmetric correlation kernel v n k n , 



Sl[p,q] = ( k pnPnP-n + kqnq n q-n + ^nfinPnq-n), 

(43) 



with the Matsubara components 



-2 yS, , ,2 

l _ \ c a°a v n , _ \ " fa^a^n 



2 2 2 

ciatvi 



(44) 



memory times r D as in Eq. (|15[) . corresponding to the 
spectral densities 



Jp,g(w) 



7p,« 1 



one simply finds 



7T 1 + T 2 U! 2 



(38) 



(39) 



in this case one exactly has k 2 = k pn k qn , and in the 
dimensionless formulation introduced in Section [II Bl thc 
coordinate pure-quantum fluctuation is given by 



a 



2 'E 



l + k v 



J (2ntn) 2 



(40) 



with the kernels given in Eq. ((18)) . For t = it becomes 
the integral 



«(0)7p)7«»t"d) 



1 + 7 P ff(ar) 



(ice 



7T x 2 + l+(7p+7 9 )5(x) 



, (41) 



with g{x) as in Eq. ((20)) . 

In Fig. [9] it appears indeed that the mean-square fluc- 
tuation of the coordinate is enhanced by the coupling 
mechanism considered in this result of the 'diago- 

nal squeezing' discussed after Eq. (|37|) . Of course, differ- 
ent distributions of the coefficients c a and d a can yield 
a less dramatic increase with the environmental coupling 
strength. 



E c a UJ a 
V 2 + UJ 2 



(45) 



Using the result of Appendix IB1 one has the expressions 
for the partition function and the mean-square fluctua- 
tions: 



1 00 2 

z = — TT ^ 

j3uj 11 v 2 (l+k n ) 2 + {a 2 +k pn 



)(b 2 +k qn ) 



(46) 



9 °° h 2 4- h 



/3 ^ ^(1+k„) 2 + (a 2 +fc p „)( fo2 + fc 9 n) ' 



<<a Q = 



Ed 2 + /Cp n 
„, x ^ 2 (l+^) 2 +(a 2 +fc p „)(& 2 +^ 



pq p 



,(47) 



while (pg + qp) = 0. Therefore, the kernel k n has the 
effect to lower the fluctuations of both canonical variables 
with respect to the case of two independent baths. It is 
reasonable to use a phenomenological guess for k qn and 
kpn, while k n is connected with them: in terms of the 
spectral densities it can be represented as 

R = 2 



LUyJJ q (u) J p (lo) 



W + vi 



(48) 



Taking Drude kernels with equal memory times r D as in 
Eq. (fl~5)) , corresponding to the spectral densities (|38|) , one 
obtains 

tin = /( r D^n) . f( x ) = ~ "T~ -1 • (49) 

Note that k n is positive and that v n k n ~ (lnn)/n — >• for 
n — > 00. In the dimensionless formulation of Section QTbJ 
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FIG. 10: System coupled with the bath coordinates and mo- 
menta. Mean-square fluctuation ol the coordinate, Eq. (|50[1 . 
at t = along the line *y P = cy q as a function of 7 9 , for differ- 
ent values of c. The dashed lines correspond to the curves for 
the two-bath case of Fig. [5] (b). 



with the kernels (I18p . the pure-quantum fluctuation of 
the coordinate reads 



a 



2 *E 



1 + K, 



^ (27rfn) 2 (l+K„)2 + (l+k pn ){l+k qn ) 



(50) 



with k n — ^/7 q 7p f(2irtnT D ). At zero temperature this 
becomes the integral 



a(0,7 P ,7?,T D ) 
dx 



(51) 



„ 7T x 2 [1+~k(x)} 2 + [l+ lp g{x)][l+j q g(x)} 



with k(x) = a/7?7p /(t d x) and g(x) as in Eq. (|20j) 

From Fig. [TO] one can see that an environmental cou- 
pling of the kind of Eq. (|42]l produces much smaller effects 
than the previously considered ones. Note that the effect 
of increasing the bath-momentum coupling j p even ini- 
tially shrinks the coordinate fluctuations, with a further 
slow rise. Moreover, it is particularly surprising to see 
that the existence of environmental coupling is irrelevant 
in the symmetric case; this means that the integral (|51[) 
sticks to the isolated system's value 1/2, whatever the 
values of 7 P = 7 g and r D : this result is far from being 
apparent from Eq. (I5T 



and momentum generally enhances their fluctuations. In 
the case of a spin system, such larger quantum fluctua- 
tions could open the possibility of reaching the quantum 
phase transition of the spin-i two-dimensional Heisen- 
berg antiferromagnet (23[, i.e., spin- lattice coupling could 
yield a vanishing zero-T sublattice magnetization with a 
power-law divergent correlation length. This possibility 
is to be further analyzed by a quantitatively accurate 
approach. Moreover and surprisingly, it turns out that 
the effect is nonmonotonic with the coupling strength, 
opening the possibility of a reentrant behavior in such 
a quantum critical system driven by its environmental 
coupling strength. 

Other kinds of coupling involving both canonical vari- 
ables and a single bath have been considered and briefly 
discussed. If the coupling terms do not simultaneously 
involve noncommuting bath variables (e.g., only bath co- 
ordinates) the result is a modification of the two-bath sce- 
nario with a further enhancement of the fluctuations and 
disappearance of the above mentioned nonmonotonicity. 
In the opposite case it clearly appears that the further 
'interference' effect gives rise to generally smaller fluctu- 
ations, i.e., the intensity of bath- momentum coupling is 
much less effective in increasing the coordinate fluctua- 
tions, which moreover show two exotic features: a coun- 
terintuitive initial decrease when 7 P is switched on and 
the intriguing stability to the isolated system's value in 
the symmetric case 7 P = 7 9 ; in other words, a symmetric 
environmental interaction appears to be ineffective for 
what attains mean square fluctuations. 
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Appendix A: The standard system-plus-reservoir 
model 



A bath attached to the coordinate q of the system of 
interest, e.g., the oscillator ([1]), is a collection of harmonic 
oscillators described by the Hamiltonian 



Hi 



+ b a (q a - q) 



^21 



(Al) 



IV. CONCLUSIONS 

Motivated by quest for a suitable description of spin- 
lattice coupling, with the lattice regarded as a harmonic 
bath of thermalized degrees of freedom, we considered 
a generalized system-plus-reservoir (SPR) model where 
both coordinates and momenta are subject to the inter- 
action with the environment. While it was known that 
coupling a bath to a canonical variable quenches its fluc- 
tuations, attaching independent baths to both coordinate 



Its characteristics are 'microscopically' specified by the 
collection of positive parameters {a Q ,6 a }, the corre- 
sponding frequencies being u a = a a b a . At variance with 
the common use [l[ we do not include here a third set of 
parameters c a in front of g, because they can be absorbed 
by means of a canonical transformation which affects the 
bath variables only, q a — > c a q a , p a — > c~ 1 p a . 

For the quantum thermodynamics at the equilibrium 
temperature T = 1//3, the influence action Si[q] corre- 
sponding to the Hamiltonian (IA1|) is obtained by tracing 
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out the bath variables, i.e., by integrating over them in 
the corresponding path integral. The influence action 
is bi-local in imaginary time (a local contribution would 
arise from a Hamiltonian operator for the system and 
couldn't describe dissipation), but it is local (or diago- 
nal) in Matsubara space [l| : 



1 



Sj[q] = / dudu 1 k{u—u') q(u)q{u') 



^ ^ fcn QnQ—n 5 



(A2) 



with the components q n for each Matsubara frequency 
v n = 2wn//3 defined by q{u) = In e- w " u . The Mat- 
subara components of the kernel k{u) are related to the 
bath parameters through 



rP u2 
k n = / du e^ u k(u) = v 2 n Y 



(A3) 



The above known results [lj qualitatively mean that stan- 
dard dissipation decreases the quantum fluctuations of 
the coordinate, while it increases those of the momen- 
tum, compared to the limit of no dissipation (k n — > 
0), where the usual quantum expressions are recovered, 
i.e., Z = psinhG&V 2 )] -1 an d (f)/b 2 = (q 2 )/a 2 = 
(l/2wo) coth(/3u;o/2). Note that in the Ohmic case, where 
k n ~ n for large n, the fluctuation of p diverges, while 
that of q remains finite. 

Of course, the symmetry of interchanging the canonical 
variables tells us that when the bath is coupled with the 
momentum, i.e., if 



(A9) 



the exact converse occurs [13H15I ]: for the harmonic os- 
cillator this is immediately clear from the fact that the 
canonical transformation p — > —q and q — > p maps the 
two coupling models onto each other. 



Starting from the Hamiltonian % = Us + Hi, after elim- 
inating the bath variables from the equations of mo- 
tion one finds that the dynamics is ruled by a quan- 
tum Langevin equation Q , where dissipation is described 
in terms of a (retarded) memory function "/(t) whose 
Laplace transform is expressed in terms of the micro- 
scopic bath parameters as 



f°° h 2 
7 (z)= / dt e -^(t) = z J2^f~2 

JO „ Z T W Q 



(A4) 



hence, it is apparent that the kernel is related to the 
memory function, 



k n = \v n \j(z=\v n \) 



(A5) 



This equality is a valuable relation, as it connects the 
quantum dissipative kernel with a phenomenological 
quantity that can be experimentally accessed in a dis- 
sipative dynamical process. 

The full system's action, i.e., the sum of the action 5s 
of the isolated oscillator and the influence action Si, is 

S = [2VnPnq-n+a 2 pnP-n + (b 2 + k n )q n q- n ], 

n 

(A6) 

and using the general result of Appendix [B] one finds the 
known results [lj for the partition function, 



-. oo 

z = — TT 



Pu ^ f 2 +a 2 (b 2 +k n ) 



(A7) 



and for the mean-square fluctuations of momentum and 
coordinate, 



V) = -rY. 



b 2 + kn 



(3 ^ u 2 +a 2 (b 2 +k n ) 



(A8) 



Appendix B: The general Gaussian integral 

Let us consider the following general quadratic action 
in Matsubara space, 

S = ~- ^2 [ a nPnP~n + 1% 9nQ-n + 2 ( 7r n+ <*n)Pn?-n] > 
n 

(Bl) 

where a n = a_„, b n — 6_„, and n n = 7r_„ arc positive, 
while the antisymmetric part S n = — <5_„ is such that 
&n ~ Vn, the Matsubara frequency, for large n. The 
Gaussian average of any quantity 0({p n , q n }) is 



(O) 



in [ ^ o({ Pn , qn} ) ^ 



(B2) 



where J = Yin^LiiP^n) 2 is the Jacobian of the Matsubara 
transformation, Z is the partition function, obtained for 
0=1, and the integrals are defined as ordinary integrals 
over the independent real and imaginary parts as 

J] f dp n = [dpo f[ 2 f dpi dp l n , (B3) 

where p± n =p^^cip n , and the analog for the coordi- 
nates. 

The action can be easily worked out by means of simple 
transformations that preserve the measure: 

(i) the coefficients of momenta and coordinates are 
made equal to u) n — a n b n by 



Pn = \ Pn , 



, n (1) 

q n = \ t~ w > 

Or, 



(B4) 



(ii) the terms that multiply 7r„ are diagonalized to 

■ji) (: 
Pn PI 



- m - (2) q^q% by the rotation 



(i) 

n 



P 



J_/l -1 
M 1 1 



Pn 



(B5) 
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(Hi) so that a further balancing is in order, 



Pn . 

V Ul n + 7T„ 

eventually giving 



p n , qh 



" ,/'^±^^. (B6) 



§ £ («„ p>~L„ + «^_ n ) (B7) 

n 

with Vt 2 n = ul-nl = albl-nl and p' n = p n -(6 n /Q n ) q n ; 
note that for f2„ to be real one must require 

7T„ < w n . (B8) 

It is now straightforward to obtain the variances 

1 n„ 



5 2 + fl 2 



(PnP-n) = (q n q-n) 



ps 2 n + n 2 n 



. , 1 5„ 
{p n q-n) = -X 



+ n-, 



(B9) 



Integrating over and g„ one finds the partition func- 
tion 



1 00 2 
/3^ 11 ^ + ^2 



(BIO) 



It is easy to go backwards through the transformations 
(Hi), (ii), and (i) to obtain the explicit result for the 
original Matsubara variables of Eq. (|B1[) , 



(PnP-n) 
(qnq-n) 

(Png-n) = 75 w I "^7" _9 ■ ( B11 ) 



p 51 + a 2 n b 2 n - 7T 2 

1 < 

P 51 + a 2 n bl - tt2 
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